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Abstract 

^ We give a detailed survey of results obtained in the most recent half decade which led to a 

O deeper understanding of the random displacement model, a model of a random Schrodinger 

^ operator which describes the quantum mechanics of an electron in a structurally disordered 

W) medium. These results started by identifying configurations which characterize minimal en- 

ergy, then led to Lifshitz tail bounds on the integrated density of states as well as a Wegner 
estimate near the spectral minimum, which ultimately resulted in a proof of spectral and 
dynamical localization at low energy for the multi-dimensional random displacement model. 
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1. Introduction 



^ By the random displacement model (RDM) we refer to a random Schrodinger operator of 

Cd the type 

B^ H^ = -A + V^, K.(x) := ^g(x-n-w„) (1) 

K^ in L'^iW^), d > 1. The potential is generated by randomly displacing translates of the single- 

\0 site potential q from the lattice sites n G Z'^. More detailed assumptions on q and the random 

9^ displacements Un will be introduced below as needed. 

(3) The RDM has proven to be much harder to analyze mathematically than the (continuum) 

1^ Anderson model 
O i7,^(^) = -A + E ^«^(^ - ^) (2) 

^ with random coupling constants A„ = A„(u;). A fundamental technical difference between 

the RDM and the Anderson model lies in their monotonicity properties. If the single-site 
potential q is sign-definite, then the Anderson model is monotone in the random variables A„ 
in quadratic form sense. This is not true for the RDM, independent of sign-assumptions on q. 
Many of the rigorous tools which have been developed to study the Anderson model rely on 
its monotonicity properties. In particular, this is true for most of the proofs of localization for 
the Anderson model near the bottom of its spectrum. In fact, if one considers the Anderson 
model with sign-indefinite single-site potential q, and thus looses monotonicity, then localiza- 
tion results are much more recent and far less complete than for the case of sign-definite q, 
e.g. [40, 32, 24, 34, 35]. The difficulties which arise are in many ways similar to the problems 
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encountered in the RDM. Related phenomena and difficulties also arise in discrete alloy-type 
models with sign-indefinite single site potential, as recently reviewed in [14]. 

Among models for continuum random Schrodinger operators, the structural disorder de- 
scribed by the RDM can be considered as physically equally natural as the coupling constant 
disorder in the Anderson model. Another natural model for structural disorder is the Poisson 
model 



^ + E 



q(x - A',), (3) 



with Xi denoting the points of a (i-dimensional Poisson process. The RDM as well as the 
Poisson model were introduced early on in the mathematical literature on continuum random 
Schrodinger operators, e.g. [26, 36] and references therein. However, progress has been much 
more limited than for the Anderson model due to the technical difficulties which arise. 

An exception is the case d = 1, where localization throughout the entire spectrum has been 
proven for the RDM and the Poisson model in [39, 8, 13]. This was possible based on the 
powerful dynamical systems methods available to study one-dimensional random operators, in 
particular those allowing to prove positivity of Lyapunov exponents and to deduce localization 
from this. However, the non-monotonicity of the RDM and the Poisson model has visible 
consequences already in the one-dimensional case, for example through the appearance of 
critical energies in the spectrum at which the Lyapunov exponent vanishes and, in some 
cases, weaker results (e.g. on dynamical localization, which has not been shown for the one- 
dimensional Poisson model). 

In dimension d > 2 it is generally expected that "typical" random Schrodinger operators 
have a localized region at the bottom of the spectrum, at least if the latter corresponds to 
a fluctuation boundary of the spectrum, which describes a boundary characterized by rare 
events. The history of localization proofs for the multi-dimensional RDM and Poisson model 
is told very quickly. For the Poisson model in d > 2 localization at the bottom of the spectrum 
was finally proven in [18] for positive single-site potentials and in [19] for negative single-site 
potentials. In both cases the powerful extension of multi-scale analysis developed by Bourgain 
and Kenig in [6] was used as a tool. 

There were two previous results on localization for the multi-dimensional RDM, [31] and 
[22]. In [31] a semiclassical version of (1) is considered and localization near the bottom of 
the spectrum is established for sufficiently small values of a semiclassical coupling parameter 
at the Laplacian. [22] considers the RDM with an additional periodic term V^gr and estab- 
lishes localization for generic (but non-zero) choices of Vper- In both works the values of the 
displacements Un have to be sufficiently small and first order perturbation effects (such as 
a monotonicity property of Floquet eigenvalues of —A + Vper in [22]) are exploited. What 
makes the "naked RDM" (1) more difficult to handle is that, as will be pointed out below, 
one ultimately has to resort to second-order perturbation effects. 

The goal of this work is to give a detailed survey of new results for the RDM (1) obtained 
in the papers [4, 5, 35] and [30], which allowed to understand that the spectral minimum of 
the RDM is a fluctuation boundary under a natural set of assumptions not requiring addi- 
tional parameters or smallness of the displacements (other than a non-overlap condition), and 
ultimately led to a proof of localization in this setting in [30]. 

The strategy used to prove localization in these works is the one provided by the Frohlich- 
Spencer multi-scale analysis [17], as described for continuum models in very accessible form in 
the book [38], and with state-of-the-art results shown in [20] and surveyed in [29]. In essence. 
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the MSA approach shows that locahzation, spectral as well as dynamical, can be proven once 
a smallness result ( "Lifshitz tails" ) for the integrated density of states at the bottom of the 
spectrum and a Wegner estimate are available as input. 

Therefore much of our effort is aimed at proving these two ingredients. However, for the 
RDM (1) one first needs to address a preliminary problem: Which configurations u = (un) 
characterize the minimum of the almost sure spectrum of H^^l To explain that this is a 
non-trivial issue, let us compare with the Anderson and Poisson models. In the Anderson 
model (2), due to monotonicity, the spectrum is minimized by choosing all coupling constants 
\n minimal (in the support of their distribution) if q is positive, while all A„ should be 
chosen maximal if q is negative. For the Poisson model (3) the spectral minimum is if g is 
positive, corresponding to regions with widely separated Poisson points. If q is negative, then 
regions of densely clustered Poisson points lead to spectral minimum — cxd. The mechanism for 
generating the spectral minimum in the RDM is much less apparent (with similar difficulties 
arising for the Anderson model with non sign-definite single-site potential). In fact, while 
for the (definite) Anderson and Poisson model the spectrum is minimized by minimizing the 
potential, for the RDM we will see that a much more subtle interaction between kinetic and 
potential energy determines the spectral minimum. 

In terms of assumptions to be made, the most important one is that the single-site potential 
shares the symmetries of the underlying lattice, here Z°'. It is fair to say that in our approach 
symmetry replaces the lack of any apparent monotonicity properties of the model, ultimately 
allowing to identify more delicate monotonicity properties which are at the core of our proofs 
of Lifshitz tail bounds and a Wegner estimate for the RDM. 

We find it remarkable how many mathematical ideas and tools had to be invoked and how 
all this ultimately fit together quite perfectly to lead to a localization proof for the RDM 
(1). Getting this across to the reader is our main motivation for providing this expository 
account of our work. Beyond merely stating a series of results, we include frequent discussions 
of the underlying motivations, often going beyond what we have been able to include in our 
previously published work. We have also tried to include at least outlines of all proofs, even 
if we frequently have to refer to the original papers for additional details. 

A rough outline of the contents of the remaining sections of this paper is as follows: In 
Section 2 we reveal how the spectral minimum of the RDM is found. In Section 3 we show 
how the proof of this is reduced to a spectral minimization property of a related single-site 
Neumann operator. This operator and its ground state properties are central to almost all our 
results. In particular, we will revisit this operator in Section 6 and explain why we ultimately 
needed to know more about it than what is stated in Section 3. To avoid having to interrupt 
the telling of our localization story, we outline the proofs of these results in Section 10 near 
the end of the paper. 

The rest of the localization story is told in Sections 4, 5, 7, 8 and 9. Section 4 yields 
information on uniqueness of configurations characterizing the spectral minimum which is 
necessary for the proof of Lifshitz tail bounds in Sections 5 and 7. The results in the latter 
two sections work in form of a boot-strap, starting with a Lifshitz tail bound under strong 
additional assumptions which are then relaxed. Our Wegner estimate for the RDM is presented 
in Section 8. In Section 9 we state the exact form of our result on localization for the RDM 
and provide references to the literature on multi-scale analysis, which show how this is proven 
based on the Lifshitz tail and Wegner bounds. In the very last Section 11 we discuss some 
open problems related to our work. 
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2. The Spectral Minimum of the RDM 

We will always assume that the displacement parameters uj = (wn)neZ'^ ^ire independent, 
identically distributed M'^-valued random variables. Their common distribution is a Borel 
probability measure on M*^. As usual, we define its support by 

supp/i := {a e M'^ : iJ,{B^{a)) > for all e > 0}, 

which is a closed set. The i.i.d. random variables a;„ can be realized as the canonical projections 

a; I— )■ a;„ in the infinite product probability space 

Under weak assumptions on fi and the single-site potential q the RDM Hi^ is self-adjoint on 
the second order Sobolev space in L'^{M.'^) and ergodic with respect to shifts in Z'' in the sense 
of, e.g., [10]. Thus its spectrum is almost surely deterministic: There exists S C M such that 

a{H^) = S for P-almost every u. (4) 

In fact, one has 

S= U aiH^), (5) 

where 

Cper '■= {w G Z'^ — >■ supp /i periodic with respect to a sub-lattice of Z'^}. (6) 

This follows by the same methods which have been used to prove a corresponding result for 
the Anderson model: That S is contained in the right hand side of (5) follows by approxi- 
mating any given random configuration with periodic configurations, truncating the random 
configuration to large cubes and periodically extending from there. On the other hand, one 
can show that almost every random configuration comes arbitrarily close to any given periodic 
configuration on arbitrarily large cubes, which is he idea behind the reverse inclusion. For a 
detailed proof, written for the case of the Anderson model, we refer to [26]. 
In particular, (5) implies that 

Eq := minS = inf{min a{Hi_^) : u G Cper}- (7) 

It is a non-trivial question to decide if there is a periodic minimizer, i.e. if the infimum in 
(7) is a minimum. In fact, we do not believe that this is true in general. Our choice of the 
following assumptions on q and fi is mostly motivated by the fact that they allow to find a 
periodic minimizer for (7). 

(Al) The single-site potential g : M'^ — )■ M is bounded, measurable and reflection-symmetric 
in each variable. Moreover, suppg C [— r, r]'' for some r < 1/2. 

(A2) Let dmax '■= k~''^ ^^^ ^ '■— {{^dmax, ■ ■ ■ , ^dmax)} dcuotc the 2'^ corners of the closure 



G oiG := {-dmax, dmax) ■ Then 



C C supp/U c G. 



The two support assumptions on q and /x have a simple geometric interpretation for the 
RDM (1): The support of each single-site term g(- — n — ujn) stays in the unit cell centered at 
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n, while it is allowed to "touch" the boundary of the cell. In fact, with positive probability 
the single-site potentials may move arbitrarily close to each corner of their cell. For a typical 
configuration of the Un see Figure 1 (where the support of q is drawn radially symmetric for 
aesthetic reasons). 



^ O ^ o 



Figure 1. The support of K; for a typical u. 

We can now identify a periodic minimizer for (7), stating a result from [4]: 
Theorem 2.1. Assume (Al) and (A2) and let u* = (ci;*)„g2d be given by 



OJ^ 



{{-lYHma.,...,{-lT'dma.), n = (m, . . . , TZrf) G Z'^. (8) 

Then Eq = rain a{H^*). 

The potential K;* (x) = J2n l{^~^~'^n) is 2-periodic in each direction and locally consists of 
densest clusters of 2'^ single-site terms placed into adjacent corners of their cells, see Figure 2. 







Figure 2. Support of K;* for d = 2. 

One may understand this result heuristically by the following strategy to construct test- 
functions which minimize the quadratic form of Hi^, at least if the single-site potential q is 
negative: The clusters in V^* form wide wells. In these wells one can place localized test 
functions with relatively small derivative, due to the width of the wells, i.e. small cost in 
kinetic energy. This gives lower total energy than the narrower wells given by individual, 
spatially separated single-site potentials. This is not how Theorem 2.1 is proven. We should 
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also point out that Theorem 2.1 does not impose any sign-restrictions on g, and thus can not 
be fully explained by the above heuristics. But the heuristics make clear that the spectral 
minimum of the RDM is determined by a non-trivial interplay between kinetic and potential 
energy. 

Instead we will give a proof of Theorem 2.1 at the end of the next section, based on the 
answer to a minimization problem for a single-site Neumann operator associated with the 
RDM. 



3. The Neumann Problem 

Theorem 2.1 provides the answer to an optimization problem involving infinitely many 
displacement parameters a;„, n e 7/. However, due to the symmetry assumptions on g, it 
turns out that the proof can be reduced to a related problem involving the optimal placement 
of just one single-site term. 

For this purpose, let Ai := {—\, \Y be the unit cube centered at the origin and — A^ the 
Neumann-Laplacian on L^(Ai), i.e. the unique self-adjoint operator whose quadratic form is 
/a |V/(a;)p(ia; for / G H^{Ki), the first order Sobolev space. 

For q as in (Al) and a G G let 

and Eo{a) := mma{H^ (a)) the non-degenerate lowest eigenvalue of Hj^ (a). For a general 
discussion of properties of operators of this type see Section 2 of [4] . 

We ask for the optimal placement of a G G to minimize EQ{a) and arrive at the following 
result. 

Theorem 3.1. Under assumption (Al) one of the following two alternatives holds: 
(i) EqIo) is strictly maximized at a = and strictly minimized at the corners C of G. 
(a) Eo{a) is identically zero. 

A proof of Theorem 3.1 under the given assumptions can be found in [4]. We will not 
discuss details of this proof here, as we will later need a strengthened version of Theorem 3.1 
for which we will also use somewhat stronger assumptions, see Theorem 6.1 and assumption 
(Al)' in Section 6 below. 

The proof of Theorem 3.1 in [4] shows that in each of the two alternatives more can be said: 

In case of alternative (i), the function E^lai, . . . ,aj, . . . , a^) is symmetric and strictly uni- 
modal in each variable. Thus, with all other variables fixed, for each j it is a strictly increasing 
function of Oj in [—dmax,0] and strictly decreasing in [0,c?max]- 

On the other hand, if alternative (ii) holds, then the strictly positive ground state eigen- 
function Uo{x,a) corresponding to Eo{a) is constant near the boundary of Ai (and thus, by 
analyticity, constant in the entire connected component of Ai \ suppg(- — a) containing the 
boundary of Ai). This reveals a mechanism which can be used to construct non-trivial exam- 
ples (with non-vanishing q) where alternative (ii) happens: 

Let 4>{x) be a positive sufficiently regular function which is constant near the boundary of 
Al and then define the potential by setting 

A(f)(x — a) 

q{x -a) = -p ^ 9 

(p[x — a) 
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as long as suppg(- — a) C Ai. Then E^la) vanishes identically for a E G and (f){x — a) is 
the corresponding eigenfunction. As follows from the proof of Theorem 3.1, this is the only 
mechanism which leads to alternative (ii). 

Alternative (i) certainly happens for all non-vanishing sign-definite potentials q, as it follows 
by perturbation theory that in this case the zero ground state energy of the Neumann 
Laplacian is pushed either up or down. But alternative (i) is generic also for sign-indefinite 
potentials, as alternative (ii) will be broken by typical small perturbations of the potential. 

Among previously known results, the ones most closely related to Theorem 2.1 can be found 
in [23] which considers similar questions for the case of the Dirichlet Laplacian — A^ instead 
of — A^. Also using symmetry assumptions on q, it is found there that the optimal placement 
of the potential in — A"^ + q{x — a) depends strongly on the sign of q. For cubic domains, 
a special case of the domains considered in [23], it is found that for positive potential the 
lowest eigenvalue is minimized if the potential is placed in a corner of the cube, while negative 
potentials should be placed into the center of the cube. This distinction does not happen in 
the Neumann case, where it is generally true that "bubbles tend to the corners" . 

While not used in our proof of Theorem 3.1 or in the proofs in [23], one can understand 
this distinction by perturbative arguments. For this, consider —A -|- \q{x — a) on L^(Ai) for 
small coupling, with either Dirichlet or Neumann boundary condition. If EQ{a,X) denotes 
the smallest eigenvalue in the Dirichlet case, then by first order perturbation theory, 

9A<(a,0)= f q{x - a)\ip{x)\^ dx, (10) 

where (p{x) is the normalized ground state of the Dirichlet Laplacian on Ai, i.e. (p{x) = 
CY[j=iCos{TrXj). In the small coupling regime minimizing (10) over a indicates the optimal 
placement of the potential. If q is positive, then the bubble should be placed into a corner of 
Ai, where l^jp has the smallest mass. On the other hand, for negative q the bubble should be 
placed into the center where the mass of l^jp is largest. 

For the Neumann case the heuristics given by first order perturbation theory is inconclusive. 
The ground state of the Neumann Laplacian is constant, and thus OxEqIq, 0) is independent 
of a. 

However, one gets correct heuristics by going to second order perturbation theory. We have 
(for a derivation see Section 2.3 of [4]) 

dMia, 0) = -2 y (^o^^(--^>^)\ (11) 

^ E,-Eo 

Here = Eo < Ei < E2 < ■ . . are the eigenvalues of the Neumann Laplacian and Uk the 
corresponding eigenfunctions. In d = 2 (for simplicity) we have that the first excited state is 
twice degenerate, Ei = E2 = vr^. Considering only these two terms in (11) (the third term 
would still give the same result) we get that dfEQ^a, 0) is approximately given by 



4 



2 / r. \ 2 

q{x — a\,y — a2) sin(7ra;) dxdy \ + I / q{x — ai, y — 02) sm{ny) dx dy 



which is non-positive. If q is reflection symmetric, then both integrals are zero for a = 0, 
indicating the position with highest ground state energy in the small coupling regime. If we 
also assume that q is of fixed sign, then both integrals become maximal in absolute value if a 
is located near one of the corners of the cube. These are the positions where the ground state 
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energy of — A^ + \q{x — a), A ~ 0, is minimal. As opposed to the Dirichlet case, the answer 
suggested by second order perturbation theory is the same for positive and negative q. 

Let us finally start to get beyond heuristics and show rigorously that Theorem 3.1 implies 
Theorem 2.1: 

Proof, (of Theorem 2.1, [4]) For any given configuration u, the restriction of if^^ to the unit 
cube centered at n G Z'^ with Neumann boundary conditions is unitarily equivalent (via 
translation by n) to H^ (wn), defined as in Theorem 3.1. Thus, by Neumann bracketing and 
Theorem 3.1, 



Tn.ma{H^) > mina 



> 



inf \^Eo{a) : a E 
Eo{a*), 



UJr, 



'(^maxj ^max\ 



where a* = {dmax, ■ ■ ■ , dmax) is one of the corners C of G. This holds for arbitrary configurations 
u and thus, by (4), Eq = minE > Eq^q*). 

Now consider u* = (wi^)„6Zd cis given by (8). The corresponding potential K;*(x) = 
Xlnezd li^ — n — u;*) is 2-periodic in Xj for each j. By Floquet-Bloch theory [37] the bottom 
of the spectrum of H^* = — A + V^* is given by the smallest eigenvalue Eq"^ of its restriction 
with periodic boundary conditions, see Figure 3. 



to Ag 



1 3\d 

2' 2' 







A^ 





Figure 3. The period cell of V^,* in d = 2. 

On Aq the potential K;. is symmetric with respect to all hyperplanes Xj = 1/2, i = 1, . . . , c?. 
Thus Eq^^ coincides with the smallest eigenvalue of the Neumann problem on Aq. Again by 
symmetry of the potential, the latter coincides with the smallest eigenvalue of the Neumann 
problem on Ai. As uJq = a*, this eigenvalue is Eq^q*). Together with (7) we have shown that 

Eo<mma{H^,)=Eo{a*). 

Combined with the previous observation that EQ{a*) < Eq this shows Eq = inina{Hi^*). D 

4. Uniqueness of the Periodic Minimizer 

After resolving the preliminary problem of characterizing the spectral minimum of the RDM, 
we could now turn to the other essential ingredients into a localization proof, a Lifshitz tail 
bound on the IDS and a Wegner estimate. However, a first look at this quickly demonstrates 
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that we also need to address the question of uniqueness in Theorem 2.1. Other than translates 
of u*, are there more periodic configurations which have the same spectral minimum? 

To motivate this, let us include a first discussion of Lifshitz tails. For this one considers 
restrictions H^l of H^ to L^(Aj^), where A^, is a cube of side-length L centered at the origin. 
As boundary condition one can generally choose what is most convenient in a given model, 
for us this will be Neumann conditions. By a Lifshitz tail bound we mean a result which says 
that the probability of H^^^l to have an eigenvalue close to Eo, the minimum of the infinite 
volume spectrum, is exponentially small in L. The meaning of "close" will be made more 
precise later. 

If u coincides with u* on Al or is very close to it, this will give a low lying eigenvalue of 
Huj,L- Our chances of getting a useful Lifshitz tail bound would worsen if there are many 
other periodic configurations with the same spectral minimum as u*, as this would increase 
the probability that random configurations are close to one of the minimizing configurations 
on A^, and thus have low lying eigenvalues. 

From this it is immediately clear that for all further considerations we will have to assume 
that alternative (i) of Theorem 3.1 holds, as under alternative (ii) it follows that Hi^^l has 
spectral minimum Eq for every configuration u. The following result is taken from [5]. 

Theorem 4.1. Assume (Al), (A2), alternative (i) of Theorem 3.1, d > 2 andr < 1/4. Then 
u* as given by (8) is, up to translations, the unique periodic configuration with m.m(j{H^^*) = 

Eq. 

Two additional assumptions were made here which deserve comment: For the "radius" r of 
the single-site potential q we require r < 1/4 rather then just r < 1/2 assumed earlier. This 
is a technical assumption, which we need to apply an analyticity argument in the proof, see 
below. Our guess is that this assumption is not necessary for Theorem 4.1 to hold. 

However, Theorem 4.1 indeed only holds in the multi-dimensional case d > 2 In the case 
d = 1 there are many periodic minimizers, as also proven in [5] : 

Theorem 4.2. Assume (Al), (A2), alternative (i) of Theorem 3.1 and d = 1. Then an 
L-periodic configuration u = {cOn)nez, ^n+L = ^n for all n G 1^, satisfies m.ma{Hi^) = Eq if 
and only if 

(i) all Un are maximally displaced, i. e. Un = ^dmm for all n, 

(ii) L is even, and 

(Hi) in each period L equally many uOn are displaced to the left and to the right. 

It is easy to see that a periodic configuration uj with these properties is a minimizer. Let (po 
be the positive ground state of —d'^/dx^ + q{x — dmax) on (—1/2, 1/2) with Neumann boundary 
conditions. It can be shown that alternative (i) implies that 

/i:=<^o(-l/2)7^^o(l/2)=:fc. 

For a configuration satisfying (i), (ii) and (iii) of Theorem 4.2 the Neumann ground state over 
the period L is found by pasting together scaled copies of (po, compare Figure 4 for an example 
with L = 4. The number of steps up is equal to the number of steps down, which allows for 
periodic extension, showing that min o"(if^) = Eq. 

Lemma 4.4 below, which holds in arbitrary dimension, establishes the necessity of (i). The 
above construction of the Neumann ground state for the period L now shows that (ii) and 
(iii) must hold for the Neumann ground state to coincide with the periodic ground state. For 
a more detailed proof of Theorem 4.2 see [5]. 
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ky]i = kVk 




Figure 4. A 4-periodic miniinizer in d = 1. 



Theorem 4.2 has surprising imphcations for the integrated density of states N{E) of the 
one-dimensional RDM. The most extreme situation occurs if /i is the Bernoulh measure with 
equal weights at the endpoints of [—dmax,dmax], i-e. 



1. 1. 

/i = -Od H — o_ 

I c\ '-'"max ' c\ 



'121 



Theorem 4.3. Let H^^ he the one- dimensional RDM with distribution /i given by (12). Then 
there exists C > such that 

N{E) > ,^_,y ^,,, (13) 



{\og{E - Eo)y 



for E sufficiently close to Eq. 



While similar phenomena have been found for Schrodinger operators with almost-periodic 
potentials, this is the first known example of a random Schrodinger operator with non-Holder- 
continuous IDS. The density n{E) = N'{E) of eigenvalues near the bottom of the spectrum 
is even higher than for the one-dimensional Laplacian where the IDS has a square-root type 
singularity N{E) = CE^/"^ at i^o = 0. Thus the randomness has the effect of pulling more 
eigenvalues towards the bottom of the spectrum, rather than pushing them away from the 
bottom as in the more common fluctuation boundary regime described by Lifshitz tails. The 
reason behind (13) is Theorem 4.2 combined with the law of large numbers. For the symmetric 
Bernoulli distribution (12) it has very high probability that in a large even period L almost 
equally many Un take values dmax and —dmax, leading to a ground state energy very close to 
Eq. For a detailed proof of Theorem 4.3 see [5]. 

We now turn back to the original goal of this section, the proof of Theorem 4.1 on the 
uniqueness of the periodic minimizer in d > 2. For this we consider a conflguration u G Cp^r 
(as deflned in (6)) and let A be the corresponding rectangular period cell. We let -ff^^ ^^"^ 
-^(^A be the restriction of H, 



respectively, and Eo{H^^) 
facts that 



to -^^(A) with periodic and Neumann boundary conditions, 
and Eo{H^^) their lowest eigenvalues. In follows from general 



mma{H^) = Eo{H^,)>Eo{H^ 



N 



A) 



(14) 



We assume that mincr(/f(^) = Eq and have to show that, up to a translation, u coincides with 



u*. This is done in two steps, 
'he flrst step establishes th, 
Neumann conditions not only on A, but on every unit cell contained in A: 



The flrst step establishes that all a;„ sit in corners and that the ground state of H^j,^ satisfles 



Lemma 4.4. Let u be a periodic configuration with inma{Hi^) = Eq. Then Un & C for all 



n 



G Z . Moreover, in this case Eo{H^ 



,AJ 



Eo{H^^) and the ground state eigenfunction ip^. 
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of H^ f^ satisfies Neumann boundary conditions on the boundary of each unit cube A„ centered 
atneAnZ'^. 

The core of the proof of Lemma 4.4 is the following calculation, based on Neumann brack- 
eting and the characterization of ground state energies as minimizers of the quadratic form: 



Eo{H, 



N 



Xv iV^o.P + /a EneAnzd qix-n- UnMu 



L ivv^^ 



neAnZ'' 



neAnZ'* 






n 



UJnMu? L IV^, 



L \^^ 



/a I^^ 



/aIV^^ 



> 



E ^» 



neAnZ'* 



4> 

JaI^^ 



^, 



0; 



(15) 



By (14) and the assumption we conclude that Eq{H^ ^) = Eo{H^^) and that all inequalities 
in (15) are equalities. We also see that Un & C for all n G Z*^, because otherwise, given 
alternative (i), the last inequality in (15) would be strict. Finally we see from equality in the 
second to last inequality in (15) that '?/'a;|A„ is the ground state for the Neumann problem on 
A„, and thus satisfies Neumann conditions on A„. 

The second step of the proof of Theorem 4. 1 is to show symmetric matching of the bubbles, 
i.e. that in each pair of neighboring unit cells within A the single-site potentials are placed 
symmetrically with respect to the common boundary of the cells. For this we use the following 
general fact from [5], to where we refer for the proof: 

Lemma 4.5. Consider a connected open region D in M"^, d > 2 and a hyperplane P that 
divides this region into two nonempty subregions. Denote by a the reflection about P and 
assume that D fl (y{D) is connected. Let E' G M and, in D, let u be a solution of the equation 

-Au = Eu (16) 

which satisfies the condition |^ = on PDD. Then u can be extended to a symmetric function 
w on D U cr{D) which satisfies the equation —Au = Eu in this region. 

To finish the proof of Theorem 4.1 let us assume that a; is a periodic minimizing config- 
uration in which, by Lemma 4.4, all bubbles sit in corners, but that there is at least one 
non-matching neighboring pair of bubbles. Let us focus on c? = 2 and the situation in Figure 5 
(the general argument in [5] uses the same idea). Circumscribe squares around the supports 



o 


o 



Figure 5. A non-matching pair. 



12 KLOPP, LOSS, NAKAMURA, AND STOLZ 

of the two bubbles and remove these two squares from the union R of the two cells. Choose 
the resulting region as D in Lemma 4.5, and a as reflection at the center line. Here (and 
only here) we need the strengthened assumption r < 1/4 in Theorem 4.1 to make sure that 
D n o-{D) is connected. 

Let u be the restriction of the ground state of H^ j^ to D. Then —Am = Equ on D and, by 
Lemma 4.4, u satisfies Neumann conditions on the centerline P Ci D. By Lemma 4.5 ti can be 
extended to a symmetric function w on DU a{D) satisfying —Aw = Eqw. But D U cr{D) = R 
and thus w is the ground state of the Neumann Laplacian on R, implying Eq = and that w 
is constant, a contradiction to alternative (i). This completes the proof of Theorem 4.1. 

To end this section let us remark that the fundamental difference in the one-dimensional 
and multi-dimensional case lies in the possibility to smoothly match Neumann ground states 
on different unit cells. In d = 1 this can always be done by re-scaling as explained by Figure 5. 
In higher dimension the boundary of cells has much more structure (is not just a point), which 
ultimately results in matching of ground states only being possible in the trivial reflection- 
symmetric case. 

5. Special Lifshitz tails 

After clarifying uniqueness questions in the previous section we finally have enough back- 
ground information to enter into a discussion of Lifshitz tail properties of the IDS for the 
random displacement model. In order to use our earlier results we will have to assume from 
here on that alternative (i) of Theorem 3.1 holds and that d > 2. We strengthen (A2) to 
require r < 1/4 and will in this section also assume that q is continuous to make use of results 
in [35]. 

Let H^]^ be the restriction of H^ to A^ = (— L — 1/2, L + 1/2)^^ with Neumann boundary 
conditions. The crucial fact required in localization proofs and also in the proof of Lifshitz 
tail asymptotics of the IDS is that the probability of H^j^ having an eigenvalue close to Eq is 
very small. The first proof of this was given in [35], where it follows as a special case of a more 
general result. The methods derived in [35] require to assume in (A2) that supp/i is finite. 
The methods developed later in [30], which are described in Section 7 below, have allowed to 
remove this additional assumption on /x. However, Theorem 5.1 is crucial as it will serve as 
the anchor for a bootstrap argument in Section 7. For this it will be sufficient to start with 
the case supp fi = C, i.e. a displacement model where all bubbles sit in corners. 

Theorem 5.1. Let the assumptions listed at the beginning of this section be satisfied and also 
assume that supp ^ = C. Then there exist C > and /i > 1 such that for all L G N, 

P (min a«J < Eo + §) < (21 + l)'^- V^'^ (17) 

In the remainder of this section we will discuss the argument from [35] which proves The- 
orem 5.1, taking some advantage in presentation from only looking at the specific situation 
which is of interest to us here. However, we will refer to [35] for many of the core analytical 
parts of the proof. 

The argument starts with decomposing the cubes A^ into quasi-one-dimensional tubes (see 
Figure 6) and restricting H^j^ to these tubes under insertion of additional Neumann boundary 
conditions. Thus let 

Al := {p e Z'^-^ : -L < pj < L foi j = 1, ... ,d - 1} 
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Figure 6. A quasi-one-dimensional tube 



and, for each p G A/,, 

L 

Ei,p:= U Ai((p,fc)), 



with Ai((p, k)) denoting unit cubes centered at (p, A;) G Z . By H^j^^ we denote the restriction 
of Hi^ to LF'iT.L^p) with Neumann boundary conditions. By Neumann bracketing we have 

peAL 
and therefore 

mina«J > min mina^^,,). (18) 

P(^Al ' '^ 

For a given u such that Ci;„ G C for all n we will say that two neighboring unit cubes are 
matching if the single site potentials in the two cubes are mirror images under reflection at 
the common boundary of the cubes. For each p G Al consider the event 

^L,p := {uj : T,L,p contains at least one neighboring pair of non-matching cubes} 

and let Xl := HpfzA^X^p. 
By independence we have 



{u i Xl,p) < fi' 



-2L 



where fi := 1/ max^ec /^(o-) (after cOp-L is chosen, the other 2L values of Up^k are determined 
by matching). This implies that 

HXl) > 1 - (2L + l)'^- V"^^- (19) 

If w G Xl, then each tube ^l,p contains at least one non-matching pair of neighboring 
cubes. Thus by Proposition 5.2 below min a{H^^) > Eq + C/L'^ for a C > independent of 
p and L. Thus Theorem 5.1 follows from (18) and (19). 

In the proposition which was used here we can without loss consider p = 0: 
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Proposition 5.2. There is a constant C > 0, independent of L, such that for every u with 
ujn (z C for all n and at least one non-matching pair of cubes in S^^q ^^ holds that 

mma{H^,L,o)>Eo + ^. (20) 

This, at least essentially, is a result proven in [35] . We will not reproduce the details of the 
proof, which involve a surprisingly rich combination of analysis tools such as a Poincare-type 
inequality, the so-called ground state transform and properties of a Dirichlet-to-Neumann 
operator, as well as some combinatorics. But we will outline the main idea: 

The proof of Proposition 5.2 starts by first arguing that is suffices to assume that 

(1) the first two cubes Ai((0, —L)) and Ai((0, —L + 1)) are non-matched, while 

(2) all other neighboring pairs are matched. 

Seeing this is not entirely trivial and requires a trick as well as some combinatorics. The 
trick consists in extending the operator H^j^ q by refiection to a twice longer tube and to 
consider the resulting operator as an operator on the torus [—1/2, 1/2]'^"^ x (M/2(2L + 1)Z) 
with periodic boundary conditions. Due to symmetry this operator has the same spectral 
minimum as H^j^ q. Now one argues that the torus can be decomposed into subsegments each 
of which has a non-matching pair of cubes at one end and otherwise only matching pairs of 
cubes. Introducing additional Neumann conditions on the subsegments lowers the spectrum 
and it now suffices to prove the claim for each subsegment (which has length bounded by 
2(2L -l- 1)). Justifying that this decomposition is possible "is an easy combinatorics, though 
somewhat lengthy to write down using symbols", where we use the words of [35] and omit the 
details. 

Under the additional assumptions (1), (2), Proposition 5.2 is essentially a special case of 
Theorem 2.1 in [35]. While our situation does not satisfy the exact symmetry assumptions 
of Theorem 2.1 in [35], the construction of our potential via matching cubes allows to mimic 
the proof in [35] almost line by line. (We mention, however, that Section 4 of [35] provides a 
slightly different argument which allows to directly use their Theorem 2.1 to prove Theorem 5.1 
above.) 

Here the results of Section 4 enter as follows: As the first pair of cubes Ai((0,— L)) and 
Ai((0, —L + 1)) does not match, the lowest eigenvalue of the Neumann problem on the union 
A2 of these two cubes is strictly larger than Eq. This follows from the arguments in the proof 
of Theorem 4.1, more specifically the argument at the end of Section 4 which ruled out that in 
this situation the lowest eigenvalue can be equal to Eq. It is this operator which plays the role 
of the operator P^ in [35]. By the results of Section 4 it is clear that inserting an additional 
Neumann condition along the surface separating A2 and S^^ \ A2 strictly lowers the ground 
state energy of H^j^q to Eq. The meaning of Theorem 2.1 of [35] is that it provides the 
quantitative lower bound CL~^ on how much the energy is lowered depending on the length 
of the attached tube. 

The value of the constant C > in (20) which is provided by the argument in [35] will depend 
on the choice of wo,-l and ujq-l+i-, after which the remaining values of Un are determined by 
matching. However, as only the finitely many values in the corners C are allowed, one can 
ultimately choose the smallest of finitely many values of C . 

As indicated earlier. Theorem 5.1 is the result which will be used later, rather than its 
consequences for the IDS N{E) of H^. However, we mention that the following Lifshitz tail 
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bound can be derived from Theorem 5.1 with standard arguments, see [35]: 

,. log|logiV(E)| ^ 1 

hmsup— — -— zr-— < — . 21 

The Lifshitz exponent —1/2 obtained here is hkely not optimal. One would expect that the 
correct exponent is —d/2, as known for the Anderson model or Poisson model. The reason 
for the discrepancy lies in the essentially one-dimensional argument which enters the proof 
through the decomposition of cubes into quasi-one-dimensional tubes. 

We stress the fundamentally different low energy behavior of the IDS of the RDM in the 
one-dimensional and multi-dimensional settings. If supp fi = C with equal probability for all 
corners, then (21) shows that the IDS has a very thin tail for d > 2, while by Theorem 4.3 it has 
a very fat tail (and thus the spectrum does not have a fluctuation boundary) for d = 1. From 
this point of view it is a fortunate coincidence that our main goal here is to prove localization 
in d > 2 and that, as discussed in the Introduction, localization for the one- dimensional case 
was already settled by very different methods. 

6. The Missing Link 

It is tempting to believe at this point of our work that we are halfway done with verifying 
the necessary ingredients for a multiscale analysis proof of localization for the RDM. Under 
suitable assumptions we have shown the Lifshitz-tail bound (17), so it remains to establish a 
Wegner estimate. Unfortunately, the assumption that supp/x be flnite in Theorem 5.1 makes 
us face a dilemma: Most known proofs of Wegner-type estimates, with the exception of some 
results in d = 1 [9, 13], require some smoothness or at least continuity of the distribution of 
the random parameters, due to the use of averaging techniques involving only flnitely many 
random parameters. For the multi-dimensional continuum Anderson model with Bernoulli 
distributed random coupling constants a localization proof near the bottom of the spectrum 
was enabled only recently by the powerful extension of multiscale analysis presented in [6], see 
also [2] for an extension to the case of arbitrary single site distributions and [21] for a detailed 
elaboration of the intricate ideas behind [6]. One of the main features of this approach is 
that the Wegner estimate is not established as an a-priori- ingredient, but its proof is part of 
the multiscale iteration procedure leading to localization. We mention that due to the use of 
unique continuation arguments this approach does not work on the lattice, leaving the proof 
of localization for the mult i- dimensional discrete Bernoulli- Anderson model an open problem. 

Thus, if we want to complete a localization proof for the RDM based on "traditional" 
multiscale analysis, the proof of a Wegner estimate will likely require a sufficient amount 
of regularity of the distribution fi of the displacement parameters Un- But this means that 
we also need to extend the Lifshitz-tail bound to more general distributions fi. The proof 
discussed in Section 5 above does not extend to this case, as it would require to take an 
inflmum over inflnitely many positive constants C with insufficient quantitative information 
available to guarantee that the resulting constant is strictly positive. 

It turned out that the missing link which allowed to overcome both remaining problems, 
the extension of the Lifshitz-tail bound to a larger class of distributions and the proof of a 
Wegner estimate for this class, is provided by an inconspicuous but crucial improvement on 
how bubbles tend to the corners, meaning Theorem 3.1. There it was shown that the function 
Eo{a), as long as it does not vanish identically, is strictly decreasing in each of its variables 
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away from the origin. The crucial improvement is that this decrease arises in the form of 
non-vanishing derivative. 

This is the first instance where we will have to require some smoothness of q, as differen- 
tiability of Eq{-) requires differentiability of q via perturbation theory, see Section 2.1 of [30]. 
For convenience, we will assume that q is C°°, even if much less is needed below and in the 
rest of this paper. 

As this is the last time we add assumptions on q, let us restate the full set: 

(Al)' The single-site potential g : M"^ — )■ M is infinitely differentiable, reflection-symmetric 
in each variable and such suppg C [— r, r]'^ for some r < 1/4. Also assume that Eo{a) = 
mma{Hj^ (a)) does not vanish identically in a G G. 

We now get 

Theorem 6.1. Assume (Al)'. Then for all a = (oi, . . . , a^) G G and alii = 1, ... ,d we have 

f < 0, ifai> 0, 

diEo{a) < =0, ifai = 0, 

[ > 0, if ai < 0. 

The proof of Theorem 6.1 is far from obvious (at least to us) and is best discussed in the 
larger context of considering similar questions for more general domains G. In order to not 
interrupt the presentation of our main story, i.e. the proof of localization for the random 
displacement model, we postpone this discussion to Section 10 below. But let us point out 
that the step from Theorem 3.1 to Theorem 6.1 turned out to be far from straightforward. The 
"smooth methods" behind the proof of Theorem 6.1 are very different from the symmetry- 
based operator theoretic methods used to prove Theorem 3.1 in [4] and, in particular, explicitly 
use second order perturbation theory. 

7. General Lifshitz tails 

The first of two important applications of Theorem 6.1 is that it allows us to extend the 
Lifshitz tail bound found in Theorem 5.1 to general distributions fi, not requiring finiteness 
of the support. 

Theorem 7.1. Assume that q and fi satisfy (Al)' and (A2). Then there exist Gi > and 
fj, > 1 such that 

P (^min a« J < Eo + ^^ < (2L + l)'^" V-'"- (22) 

for all L G N. 

We will prove this result by comparing the quadratic form of H^^^ with the quadratic form of 
a modified displacement model where all bubbles have been moved to the closest corner within 
their cell. Thus, for a & G, let c{a) G C be the corner closest to a (if several corners are equally 

close, any of them can be chosen). For a displacement configuration u = {uJn)n&Z'' ^ G , 

define c{uj) EG by (c(a;))„ = c{un). 

From Theorem 3.1 we know that the single-site operator Hf^_^{c{a)) has lower ground state 
energy than H^ (a). Theorem 6.1 allows us to quantify this, saying that the distance of the 
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two ground state energies is proportional to \a — c(a)|. In particular, there exists C2 G (0, 00) 
such that 

Eo{a)~Eo>^D{a), (23) 

where D{a) = miucgc |o — c\. This is one of the two central ingredients in the proof of the 
following result. The other one will be Neumann bracketing. 

Proposition 7.2. There exists a constant C3 G (0, 00) such that, in the sense of quadratic 
forms, 

1 ^ 



<L-^0>^(i^eM,L-^0). (24) 



for all u E G and all L >0. 
In particular, (24) implies 



1 



mina«J - E, > -{mma{H^^^^^^ - E^). (25) 

The RDM -ffc(i.j) has i.i.d. distributed displacements supported on C and thus satisfies the 
assumptions of Theorem 5.1. Therefore, with C and /i from Theorem 5.1, 



P ( mina(iJJ);j - E, < ^^ < P (^min a(iJX),L) - i?o < ^ ) < (2L + 1)^^V"'^ 



^] < P (mina(iJ5L),J - i^o < ^ 



proving Theorem 7.1. 

Thus it remains to prove Proposition 7.2. The strategy for this is to first prove a correspond- 
ing result for the single-site operators H^ (a) and then extend this by Neumann bracketing 
to the operators H^^ ^. For the single-site operators one separately considers the cases where 
a is close to a corner or not close to a corner. 

Lemma 7.3. There exist C > and 5 > such that, if D{a) < 6, then 

Hlia) -Eo> ^(Hlic) ~Eo + \a- c\). (26) 

Lemma 7.4. Fix 6 G (0, 1). There exists Cs G (0, 00) such that, for D{a) > 6 and all c E C, 

Hlia) -Eo> ^«(c) -Eo + \a- c\). (27) 

Before discussing the proofs of the two Lemmas, let us show how we use them to prove 
Proposition 7.2. Note that, applying Lemma 7.4 with S as provided in Lemma 7.3, both 
Lemmas combined prove the L = case. To extend this to general boxes we employ an 
argument previously used in the proof of Theorem 2.1 in [34]. It is crucial here that we work 
with Neumann boundary conditions. 

For ip G H^{A2L+i), the form domain of H^j^, one has that the restriction of tp to Ai(ra) is 
in H^{Ai{n)) for each n G Ag^,^^^ := A2L+1 H Z*^. Moreover, 

where we work with the usual slightly abusive notation for quadratic forms. 
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The same argument may be applied to Hc[uj), 

"SA^L + l 

Now Proposition 7.2 follows by applying Lemmas 7.3 and 7.4 for each tt,, summing, and 
omitting the positive term Xlnd'^n ~ c(cu„)|V', V')Ai(n)- 

Before we can end this section, we still owe a discussion of the proofs of Lemmas 7.3 and 
7.4. To see the latter, note that D{a) > 6 implies Hj^_^{a) — Eq > 6/C2 by (23). Using the 
rough bound \q{x — a) — q{x — c)\ < 2||g||oo and setting C := 1 + 2C2\\q\\oo/S we get 

{C + l){H^a) - Eo) - «(c) -E,)>^- 2||g|U = ^, 



and thus 

Hlia) -Eo> ^ (<(a) " i^o + ^ 

As 6 = r-^\\ci — cl > „, ^ ^ \a — cL (27) follows with 1/Cs chosen as the smaller of 1/(C + 1) 

aiid6/{2{C+l)C2dmax^)- 

The previous argument doesn't use the full strength of (23), but only that -E'o(') is continuous 
and strictly minimized in the corners. The proof of Lemma 7.3 is more subtle and depends on 
the linear growth of -E'o(') away from the corners. To a G G pick c eC such that D{a) = \a — c\. 
By smoothness of q we have the Taylor approximation q{- — a) — q{- — c) = {c — a) ■ Vq'(- — 
c) + o{\a — c\) and thus 

Hl{c) -Eo = Hlic) ~Eo + {c-a)- Vg(- - c) + o{\a - c\). (28) 

Bounding the left hand side by (23) we get, in the sense of quadratic forms, 

H^ic) -Eo + {c-a)- Vqi- - c) > —\a - c\ + o{\a - c\). 

<-^2 

Hence, for p G (0, 1) sufficiently small and a G S'^^^ with a = c + pa E G, 

Hl{c)-E,-pa.Vq{--c)>^^. 

We apply Lemma 7.5 below with A = Hj^^{c) — Eq and B = —pa ■ Vg(- — c) to conclude that 
for Cp = max(2, 2C2/P) t G [0, 1/2] and a G S'^'^ with c + pa eG, 

Hlic) -Eo- tpa ■ Vg(- - c) > ^«(c) - Eo + t). 

From this and (28) we find for |a — c| < p/2 and t = \a — c\/p, 

H^ia) -Eo> UhI{c) -Eo + \a- c\/ p) + o{\a - c\). 

This implies (26) if 5 > is chosen sufficiently small, completing the proof of Lemma 7.3. 
We have used the following simple fact, which was previously used in a similar context in 

[34]. 
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Lemma 7.5. Let A be self-adjoint and B hounded and self-adjoint with A > and A + B > 
Co > 0, then 

A + tB > mm(-,co) -(A + t) 

for all te [0,1/2]. 

This is elementary: 

A + tB = {l-t)A + t{A + B) > -A + tco>mm{-,co){A + t). 

8. Wegner Estimate 

To describe the ideas behind the proof of a Wegner estimate, we consider if^, a suitable 
random Schrodinger operator on L'^{W^). Let L > and H^ j^ be the restriction of H^ to the 
cube Al with, say, Dirichlet boundary conditions. The boundary conditions are expected not 
to play a too important role. 

A Wegner estimate (see [41]) is an estimate on 

E(tr X[Eo-e,Eo+e]iK,L)) = E(#{eigenvalues of i/;,^ in [^o - £, ^o + e]}) (29) 

for L large, e small and a fixed energy Eq. It can also take the form of an estimate on the 
probability P{-f/^^ j^ has an eigenvalue in [Eq ~ e,EQ-\-e]} which, by Chebyshev's inequality, is 
smaller than the previous quantity. 

From their very form, it is clear that both quantities should increase with e and with L. 
The existence of an integrated density of states for H^ suggests that the optimal upper bound 
should be proportional to \Al\ ~ L'^. The optimal upper bound in e is related to the regularity 
of the integrated density of states. The best bound one may expect is of the form CeL'^. 

Let us give the heuristic underlying such a bound in the simplest case, the case when HJ^ is 
the continuous Anderson type model Hxiui) defined in (2), when q has a fixed sign, say, positive, 
is continuous and bounded, its support contains Ai and the coupling constants (A„)„g2d are 
i.i.d. and bounded. 

Let {Ej{u,L))j denote the eigenvalues of Hx(ui) l ordered increasingly. To estimate the 
quantity E(tr X[Eo-e,Eo+e]iH^^^^L)), we can write 

E(tr X[Eo-e,Eo+s]{H^,L)) = ^ i^X[Eo-e,Eo+e]{Ej{(^,L)) 

< ^P{E,(u;,L)G[Eo-£,Eo + £]}. 

where, by standard bounds on Schrodinger operators, i^N^ < L'^. 

In the case of the continuous Anderson model under the assumptions made above, for a > 0, 
the operator inequality H^._^_- — Hx(w) ~ ^ ^^^^^ ^^ ^^^^ 

Vj, E,{\{uj) + a,L)~ E,{X{io), L) > «, (30) 

where a is the vector whose entries are all a. 

Based on this and under the assumption that the distribution of the A„ has a bounded 
density one can prove that 

Vj, F{Ej{uj,L)e[Eo-e,Eo + e]}<e (31) 
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and obtains the desired bound in eL'^. The proof of (31), while essentially based on the ideas 
described above, requires additional technical work. 

For the discrete d-dimensional Anderson model, the bound 

n^T X[E,-e,E,+e]{Hl^L)) <CeL'' (32) 

essentially goes back to Wegner's original paper [41], with some technical details filled in later. 
A detailed proof, essentially following Wegner's orginal argument, can be found, for example, 
in the recent survey [28]. Obtaining the bound (32) for the continuum Anderson model was 
harder, as one can not use the same rank one perturbation methods as in the discrete case 
to control the spectral shift due to single site terms. Initially, a bound of the form CeL'^'^ 
was obtained for the continuum Anderson model in [27] (for a proof with slightly different 
methods see also [38]). For q of fixed sign, but without the assumption that the support of q 
contains Ai and at arbitrary energy, the linear in volume bound (32) was ultimately obtained 
in [11]. 

To describe how a Wegner estimate for the random displacement model considered here 
was found, let us describe a generalization of the idea outlined above which goes back to 
[31, 32, 33]. To estimate f{Ej{u, L) G [Eq — e,EQ + e]}, we study the mapping uj \-^ Ej{u, L) 
that realizes a "projection" from the parameter (probability) space onto the real axis; and we 
want to measure the size (with respect to the probability measure on the parameter space) 
of the pre-image of some interval. The idea is then to find a vector field V in the parameters 
u such that the eigenvalue Ej{u,L) moves when u moves along the fiow of the vector field. 
The fiow of V foliates the parameter space nicely and the volume we want to measure is just 




Ej(co) 



^-e %s 



Figure 7. Foliation and projection of the probability space 



the volume contained in a layer between two leaves (see Figure 7). This volume will then be 
of size the width of this layer at least when the probability measure has a regular density. So 
if one is able to do this for all the eigenvalues, one gets an estimate of the form eL'^. 

To be able to do this for all eigenvalues at a time, one may choose V so that H^ dif- 
ferentiated along V has nice properties (e.g. positivity). Let us take the simple example 
of the continuous Anderson Hamiltonian under the assumptions made above. If we take 

d 



V = div 



X{oj) 



E 

neAi 



dXn{uj)' 



then VH^^^^ 



> 



XAr 



This ensures VEj(u,L) > 1 which is 
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exactly (30). This is what is needed for a Wegner estimate for the continuous Anderson 
Hamiltonian when the single site coupling constants admit a bounded density. 

The right choice of vector field is model dependent; for the Anderson model, discrete or 
continuous, in many cases, one may use the divergence vector field as above (for example, 
see [41, 38]). Another useful vector field with respect to this problem is the generator of 
the dilations Xln -^"('^)^a„(w)- It can be used to get a Wegner estimate for the continuous 
Anderson model without sign assumptions on V [32, 24] and, in certain cases, for the random 
displacement models [22]. Other types of randomness may require different types of vector 
fields, see e.g. [31, 15]. 

For the random displacement Hamiltonian H^^, we will use the same idea and introduce a 
new vector field. The choice of this vector field is motivated by Theorems 3.1 and 6.1, which 
indicate that at least low lying eigenvalues should decrease monotonically if the single site 
potentials are moved towards a corner of their cell. 

For a function / on G we set 

{dJ)ia):=^^f\^-Vfia), 
\c[a) — a\ 

with c{a) denoting the corner closest to a as in Section 7. Thus, dc denotes the directional 
derivative in the direction of the closest corner, where points a with multiple closest corners 
will not play a role in the arguments below (starting from (34) below we introduce a cut-off 
which restricts the values of a relevant for the proof to small neighborhoods of the corners). 
By Theorem 6.1 there exist ^o > and tq > such that 

dcEo{a) < —6o for all a G Ar^ := {a E G : |c(a) — a\ < ro}, (33) 

a neighborhood of C. 

Let rj G C°^(]R) such that < ?7 < 1, rj(r) = 1 for r < Tq and r]{r) = for r > 2rQ. Using 
this function as a cut-off, we localize the vector fields associated with dc onto a neighborhood 
of the corners, defining 

{d'J){a):=v{\c{a)-a\){dJ){a). (34) 

For each n G Z*^, we write 

d'cu.^H^ = d'c^^M- -n-Un) = -v{\c{(^n) - ujn\)^^-^^ — ^ " (Vg)(- - n - w„). (35) 

\C[UJn) — ^n\ 

If ip E H^{A2L+i), the form domain of -f/^^^, then ipn '■= ^|Ai(n) ^ -f^^(Ai(n)), the form 
domain of Hn{un), and, with the usual abuse of notation for the quadratic form, 

(^,<^^)= Yl (^n,i^n(c^n)^„), (36) 



"^^Kl+i 



as well as 



E (^'C„<L^)= E (^n,5U^"K)^-)- (37) 



"e^2L+l "-^Kl+i 



Proposition 8.1. There exist 6i > and ^2 > such that 



"e^2L+i 



for all LeN, andtpe H\A2l+i) with (^, [H^^ - Eq)^) < 62 
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Thus, near the bottom of the spectrum of H^^l, the vector field — 2. ^c,w„ satisfies exactly 

the property we are looking for to proceed according to the heuristics explained above. For the 
details of the proof of Proposition 8.1 as well as for the proof of the implied Wegner estimate 
in Theorem 8.2 below we refer to [30]. 

In order to exploit Proposition 8.1 we formulate the following final set of assumptions for 
the distribution /i of the i.i.d. random displacement parameters u = {oJn)'- 

(A2)' With G and C as above, let C C supp/i C G. Also assume that there exists a 
neighborhood of C on which /i has a C^-density. 

More formally, this means that there exists e > and a C^-function p : G — ^ M, such that 
for every S C Uaeci^ '■ \x — a\ < £:} fl G we have 

Js 

This guarantees that the random variables |c(u;„) — w^l, which parametrize the layers of 
the vector field, have absolutely continuous distribution near (i.e. for values of a;„ near the 
corners). Thus we have exactly the situation required by the above heuristics, which indeed 
leads to the following Wegner estimate: 

Theorem 8.2. Assume (Al)' and (A2)'. Then there exists S > such that, for any a G 
(0, 1), there exists Cq > such that, for every interval I C [Eq, Eq + 6] and L G N, 

HtrXiiKi)) < CMrL". (39) 

We finally need to comment on the reason for the appearance of the exponent a G (0, 1) 
in (39). This is a final price we pay for the non-monotonicity (as well as non-analyticity) of 
our model in the random parameters. The latter necessitate some additional changes in the 
original strategy of Wegner. The proof of Theorem 8.2 in [30] uses an adaptation of a Wegner 
estimate proof developed in [24] to handle Anderson models with sign-indefinite single-site 
potentials. Their argument is based on L^-bounds for Krein's spectral shift function proven 
in [12]. These bounds only hold for p < oo and not for p = oo, which is the reason that we 
can't choose a = 1 in Theorem 8.2. 

9. Localization 

At this point we have essentially reached the end of the story which was to be told here. 
As explained in the introduction, with the Lifshitz tail bound of Theorem 7.1 and the Wegner 
estimate of Theorem 8.2 available, localization for the RDM can be proven via multi-scale 
analysis. The MSA method is as powerful as the details of carrying it out are intricate. 
Presenting these details is not a goal of this article. Very good introductions into the mathe- 
matics of MSA can be found in the book [38] and the surveys [29] and [28], which also provide 
extensive bibliographies. 

The main task left to us is to state state the exact result on localization for the RDM which 
was obtained in [30]. Here Xx denotes the characteristic function of a unite cube centered at 
X, Xi{H) the spectral projection onto I for the operator H, and || ■ II2 the Hilbert-Schmidt 
norm. 
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Theorem 9.1. Assume (Al)' and (A2)'. Then there exists 6 > such that iJ^ almost surely 
has pure point spectrum in [Eq, Eq + 6] with exponentially decaying eigenf unctions. 

Moreover, H^ is dynamically localized in I , in the sense that for every C < 1; there exists 
C < oo such that 

E ( sup \\x.g{H^)Xi{H^)Xy\\l\ < Ce-l^-^l' (40) 

Vl9l<l / 

for all x,y E Z"'. The supremum is taken over all Borel functions ^f : M — > C with satisfy 
\g\ <1 pointwise. 

We note that dynamical locahzation in the physical sense is covered by (40) in choosing 
g{H) = e"**'^ and taking the supremum over t G M. 

The subexponential decay in |x — y| found in (40) is the strongest type of dynamical lo- 
calization which has been obtained through MSA. This is a result of Germinet and Klein in 
[20], who used a four times bootstrapped version of the MSA argument, allowing to conclude 
strong forms of localization from rather weak forms of initial length estimates provided by the 
Lifshitz tail bound. As described in some more detail in [30], the survey paper [29] provides 
a very useful resource in explicitly singling out all the properties of a model, which go into 
the argument in [20]. In addition to the crucial Lifshitz-tail and Wegner bounds proven in 
Theorems 7.1 and 8.2, the RDM has all other required properties. 

In this context we also recommend the book [38] as a very readable account of MSA. Similar 
to [29] it clearly exhibits the properties of a model which are needed to prove localization 
via MSA. It uses a version of MSA less sophisticated than what is done in [20], essentially 
bootstrapping the MSA scheme just twice, to conclude spectral localization and a weaker form 
of dynamical localization, that is 

E I sup \\\X\Pg{H^)xi{H^)xo\\] < oo (41) 

Vl9l<l / 

for all p > in a p-dependent neighborhood I oi Eq, giving power-decay rather than subex- 
ponential decay in (40). Here |A| is the multiplication operator by the length of the variable 

x eR'^. 

10. Bubbles Tend to the Corner 

10.1. Bubbles tend to the boundary. The Lifshitz tail estimate as well as the Wegner 
estimate have as a basic input an estimate on the ground state energy of a Neumann problem 
as a function of the position of the potential. 

It is convenient to take a more general point of view for this problem in which the unit 
cell Al is replaced by a bounded domain D with smooth boundary, and to consider the non- 
degenerate lowest eigenvalue Eq^q) of the operator 

- A^ + q{x - a) . (42) 

Here A^v is the Laplace operator with Neumann boundary conditions on dD. We shall assume 
that the potential q is smooth and has compact support such that the set 

G = {aeR'^ : supp q{- - a) C D} 

is not empty. Thus, G is an open and bounded set. We shall, in addition, assume that it is 
connected. 



24 KLOPP, LOSS, NAKAMURA, AND STOLZ 

The following second order perturbation theory result sets the stage for our investigation. 
Denote hy da = w ■ V a and likewise d^ = w ■ Vx where ly is a fixed vector and the subscript 
denotes the variable in which we differentiate. 

Lemma 10.1 (Second order perturbation theory). The lowest eigenvalue as a function 
of a satisfies the equation 



B{uk,dxUoY 



dlEo - AdaEoiuo, dxUo) = 2 / V ■ idxUoVdxUo)dx - 2 V ^r^'^'^"^ . (43) 

Jd T^o ^^- ^0 



Here B{u,v) denotes the bilinear form 

B{u,v) = {u,Av) — {Au,v) , 

and Uk{x; a) is the eigenfunction associated with the eigenvalue Ek{a), k = 0,1,2, ... . 

The proof of this lemma can be found in [4], with additional modifications in [30]. By 
Gauss's theorem, the first term on the right side can be written as 



2 / dxUoN{x)-V{dxUo)dS{x) , 

JdD 

where N{x) is the outward normal at the point x G dD. The following computations should 
reveal somewhat the geometric structure of this term. For any fixed point x we can extend 
N{x) to a smooth vector field in a neighborhood of x. We write 

N{x)-V{dxUo) = ^WjNi{x)djdiUo = ^wf{x)Ni{x)djdiUo+iw-N{x))^Nj{x)Ni{x)djdiUo , 

(44) 
where 

w^{x) = w — {w ■ N{x))N{x) 

is the projection of w onto the plane tangent to dD at the point x. For the first term on the 
right side of (44) we write 

^wf{x)Ni{x)djdiUo = -^wf{x){djNi{x))diUo + ^wf{x)dj{Ni{x)diUo) 

ij id i,j 

and note that the last term vanishes. Indeed, this term is a tangential derivative of the function 
N{x) ■ Vuq which vanishes identically on dD {uq is the Neumann ground state). Also, 

Y]w^(x)(9jA^j(x))9jMo = w^lx) ■ K{x)Vuo, 

i,j 

where K{x) = {Kj^i{x)) = {djNi{x)) is the curvature matrix of dD at the point x. Thus we 
have that 

N{x) ■ V{dxUo) = -w^{x) ■ K{x)Vuo + {w ■ N{x))^Nj{x)Ni{x)djdiUo . 

Once again, since N ■ Vmq = 0, 

dxUo = w ■ VxUo = w^ (x) ■ VxUq 



and hence 
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dxUoN{x) ■ V{d^uo)dS{x) 

dD 



- 2 / w^{x) ■ V^UoW^{x) ■ K{x)VuodS{x) 



dD 



To summarize we have 

Lemma 10.2 (Second order perturbation theory). The lowest eigenvalue as a function 
of a satisfies the equation 

dlEo - 4daEo{uo,d,uo) = -2 / w^{x) ■ V,uow^{x) ■ K{x)V^uodS{x) 

JdD 



+ 2 w^{x)-V^Uo{w- N{x))y2^ji^)^ii^)9j9i'^odS{x) 

JdD ~i 



^ Ek-Eo ' ^ ' 

This formula remains correct if apphed to a rectangular parallelepiped and a derivative 
da parallel to the edges of the domain. While the curvature matrix becomes singular along 
the corners and edges of the parallelepiped one may argue that these singularities do not 
contribute to the right hand side of (45) because the derivatives of Uq vanish in the directions 
in which K is singular. A direct argument for this case is provided in [30]. In this case no 
curvature term appears as the faces of the parallelepiped are flat. Moreover, the second term 
in (45) vanishes since one of the two terms w-^{x) or w ■ N{x) always vanishes. Thus, one gets 

Corollary 10.3. // the domain is rectangular, then for all derivatives da parallel to the edges 
of the domain we have 

dlEo - 4daEo{uo, d^uo) = -2j2 ^^^''^^^'^' . (46) 

The point of this formula is that the right side has a definite sign. Corollary 10.3 will be 
the crucial input for showing that the energy minimizing position is when the potential sits 
in the corners. 

Let us expand the scope a little bit by considering smooth domains D. For such domains, 
by summing over the canonical basis of unit vectors w one obtains the formula 

Ai^o - 4(«o, Vno) ■ VEo = -2 / Vuo ■ K{x)Vu, dS - 2^ ^i^J^^^^^f^^ . (47) 

JdD t^^ Ek-Eo 

This is best seen by only using the first identity in (44) in the above argument without 
introducing vj-^, see also [4] for a direct proof. The right side of (47) has a definite sign for the 
case where the boundary has a positive curvature matrix, i.e., is a convex surface. Assuming 
that the right side of (47) is given, this equation can be considered as a second order elliptic 
equation for the eigenvalue EQ{a) and hence it is amenable to a strong minimum principle 
(see, e.g., [16] Theorem 3 on p. 349) provided we know that the right side of (47) is strictly 
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negative. It is then a consequence of the strong minimum principle that if EQ{a) attains its 
minimum over G at an interior point of G, then EQ{a) is constant throughout G. 

One case where one may try to exploit this reasoning is if the domain D is strictly convex in 
the sense that the curvature matrix K{x) is positive definite at every point of the boundary dD. 
If the eigenvalue EqIq) attains its minimum at oq G G then V-Eo(c^o) = and A£'o('^o) ^ 0- 
Hence the right side of (47) must vanish. Since K{x) is assumed to be strictly positive at all 
points of dD we have that Vmq vanishes on dD, i.e., uo{x, ag) is constant there. Using this it 
is not hard to prove the following theorem. 

Theorem 10.4 (Strong minimum principle for Eq). //£'o(oo) = i^^a&G Eo{a) for some 
ao G G, then EQ{a) is identically zero. In this case the wave function is constant in the con- 
nected component of the complement of the support of the potential that touches the boundary 
dD . 

In other words, if -Eo(ct) is not identically zero, then it attains its minimum on the boundary 
of G. 

Theorem 10.4 corresponds to Theorem 1.4 in [4] where a proof is given. Theorem 1.4 as 
stated there is slightly inaccurate since it is implicitly assumed that the complement of the 
support of the potential has only one connected component (the exterior component touching 
dD). 

For a given domain and potential it is in general not easy to verify if the lowest eigenvalue 
is independent of the position of the potential. A construction of examples with this property, 
starting with the ground state eigenfunction, was described in Section 3, see (9). But it is 
relatively easy to verify non-constancy of EQ{a) in a number of cases. One has to exhibit one 
position of the potential where the ground state energy is not zero. This is obvious if the 
potential has a fixed sign. Likewise, if the potential is not identically zero and if its average 
is less than or equal to zero, we can use the constant function as a trial function and see that 
the ground state energy is strictly negative, for if it were zero, the constant function would be 
the eigenfunction and the potential would be identically zero. 

An interesting conclusion can be drawn using Hopf 's lemma. Since the domain is smooth, 
every point satisfies an interior ball condition. Thus, by Hopf's lemma (see, e.g., [16] p. 347) 
we conclude that the derivative of Eo{a) at the point where the minimum is attained and 
normal to the boundary is strictly positive. 

While all these ideas put us on the right track, we cannot apply them directly to our 
situation. The underlying domain Ai is not strictly convex and its boundary is obviously not 
smooth. Moreover, we need the minimal configuration to be in the corners and not just on 
the boundary, and most importantly we need estimates on how the eigenvalue increases away 
from the boundary. Note that Hopf's lemma will not do, since at the corners the interior 
ball condition is not satisfied. Nevertheless, by a refined analysis we can show that these 
statements remain true for the case where the domain is a rectangular parallelepiped. 



10.2. Eq has a non-vanishing derivatives. Let us consider formula (46) for the case where 
we take the derivative in the 1-direction. We write G = I x G' with an open interval / and 
an open d — 1-dimensional rectangle G'. As we shall fix the variables 02, • ■ ■ , a^ in G' we shall 
suppress them from the notation and just write Eo{ai). Thus, the variable ai varies over the 
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open interval / which is symmetric with respect to the origin. (46) takes the form 

Ek — Eq 



Eq - 4:Eq{uo, diuo) = -2 2^ F _ -p — • (^^ 



The following dichotomy is of interest for us. 

Theorem 10.5. Assume that the right side of (48) vanishes for some ai^ G /. Then E^lai) = 
identically in I and for every ai G / the eigenfunction UQ{x,ai) is constant in the connected 
component of the complement of the support of Qa that touches the boundary of D 

Let us give a sketch of the proof, for more details see [30] . If this function vanishes for some 
value of Oi^O) then we must have that 

B{uk,diUo)=0,k = 0,1,2,... . 

In other words 

{uk, AdiUo) = (Aufc, diUo), k = 0,l,2,... . 

This can be used to show that diUQ satisfies a Neumann condition on the faces Si and Ti of 
G perpendicular to the 1-direction, see the proof of Lemma 2.4 in [30]. Thus Ux^^xi = on 
^i and Ti. Since ai^ G / and (02, . . . ,ad) E G' the potential is zero in a neighborhood of the 
boundary dD. On the faces Si and Ti the function uq satisfies the equation —A'uq = EqUq, 
where A' = A — ^. Recalling that all the first derivatives of Uq vanish on the intersections 

of the faces (e.g. on the intersections of 5*1 and Ti with other faces of G), Uq must satisfy a 
Neumann condition on the boundary of the faces Si,Ti. Since Uq is non- negative it is the 
ground state eigenfunction of —A'uq = EqUq and hence must be constant. Thus, £'0 = 
and Mo is harmonic outside the support of the potential. If we pick a point xq on Si away 
from the support of the potential (such that the potential is zero near xq), we may assume 
by the reflection principle, that the function uo{x,aifi) is harmonic in a neighborhood U of 
xq. Moreover, uo{x,aifi) is constant on U H Si and didjUo{x,ai^o) = ioi j = 1,2 . . . ,d. 
From this one can easily conclude that Uo{x, ai^) must be constant in U (see [30]) and hence 
in the component of the complement of the support of Qa that touches the boundary. This 
proves Theorem 10.5, with the ground state being obtained by shifting the potential and Uq 
simultaneously. 

Armed with this information we can prove Theorem 6.1 easily. It suffices to consider 
the 1-direction. Returning to (48) we may introduce an integrating factor F{ai) such that 
F'(ai) = —4:{uo,diUo) and rewrite this equation as 

(e^E',y (ai) = -2e^(-) V ^^^^'^^^°^' . 

By Theorem 10.5 and assumption (AI)' the right hand side is strictly negative. Since we 
assume that the potential is symmetric with respect to reflections about the 1-direction, we 
know that the eigenvalue Eo{ai) is a symmetric function, i.e., Eo{—ai) = Eo{ai) and hence 
E'q{0) = 0. Thus, for ai > 0, by integrating we obtain 

E'iai) = -2e-^('^^) H e^(") V ^^^''^Z'^\ a < , (49) 

Jo t^o ^^~ ^0 

which implies Theorem 6.1. 
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11. Open Problems 

The main reason for presenting this expository account of our results on the random dis- 
placement model is that the methods developed have led to a satisfactory understanding of 
localization for this model, where multiple pieces of a puzzle eventually fell into place to reveal 
a complete picture. Thus, to some extend, this is the end of a story. Nevertheless, various 
aspects of our work reveal natural and non-trivial questions for further work. We end our 
presentation by describing some of them. 

11.1. Bubbles tend to the boundary. Our work has led to two types of results about the 
optimal placement of a potential to minimize the lowest eigenvalue of the Neumann problem 
(42) on a given domain G. While Theorem 10.4 shows that for general convex domains the 
minimizing position is at the boundary. Theorems 3.1 and 6.1 characterize corners as the exact 
minimizing positions for the special case of a rectangular domain. It is natural to believe that 
corners are good candidates for minimizers for other polyhedral domains, in particular for 
regular n-gons in rf = 2, but our methods do not allow to prove this for any n 7^ 4. The 
second term in (45) does not vanish unless the domain is a rectangular parallelepiped! Of 
course, in this case one would assume that the potential shares the symmetries of the domain. 
Radially symmetric potentials would be the most natural candidates and it might also help 
to choose them sign-definite. Particularly interesting would be a proof of this for equilateral 
triangles, n = 3, because all the other results presented in this paper would then apply to get 
a localization proof for the RDM on triangular lattices. 

It is natural to conjecture that on smooth convex domains the minimizing position of the 
potential along the boundary should be at a point of maximal curvature, as long as the 
potential has sufficiently small support to smoothly fit into the boundary at this point. A 
case where one could hope to formulate and prove this rigorously is an elliptic domain G in 
M? with a small radially symmetric potential. 

11.2. Periodic minimizers. Formula (5), characterizing the almost sure spectrum of the 
RDM in terms of periodic displacement configurations, holds under very weak assumptions 
on n and the single-site potential q. Is it necessary to take the closure on the right hand side? 
That we were able to prove the existence of a periodic minimizer in Theorem 2.1 was due to 
a number of lucky coincidences (or of well chosen assumptions, for that matter). One might 
ask if for more general cases the existence of a periodic minimizer in (5) is the rule or more 
likely to be an exception. 

For example, what happens if we keep all the assumptions in (Al) and (A2) above, with 
the exception that suppyU may not contain all the corners C of G? Our simple method of 
constructing a minimizer by multiple refiections breaks down. It is not clear at all, and 
maybe not likely, that a periodic minimizer exists. The same problem arises if we don't 
require that the single-site potential is refiection symmetric in each coordinate. One might be 
led to believe that periodic minimizers hardly ever exist, but for the moment we do not know 
a single counterexample. 

One may also look at non-rectangular lattices, while keeping all the desired symmetry 
assumptions on q and supp/i. If, as suggested in Section 11.1, bubbles tent to the corners for 
regular triangles, then a periodic minimizer for the triangular lattice in M^ could be constructed 
by repeated refiection, as in the case of rectangles. However, the same does not work for a 
hexagonal lattice, where contradictory positions of some of the bubbles arise after just a few 
refiections. Is there nevertheless a periodic minimizer for the hexagonal RDM? 
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11.3. The fractional moments method. In our proof of localization for the RDM we 

followed the strategy provided by multi-scale analysis, which yields spectral and dynamical 
localization based on the Lifshitz tail bound and Wegner estimate provided by Theorems 7.1 
and 8.2. Another method which has provided localization proofs for multi-dimensional ran- 
dom operators is the fractional moments method (FMM) originally introduced by Aizenman 
and Molchanov in [3] to give a simple proof of localization for the multi-dimensional discrete 
Anderson model. This method has meanwhile been extended to show localization for contin- 
uum Anderson models [1, 7]. While likely not being as universally applicable as MSA (and 
quite certainly not extendable to situations as considered in [6]), an interesting feature of 
the FMM is that, in situations where it is applicable, it yields a stronger form of dynamical 
localization than what has been obtained via MSA. Instead of (40) one gets the stronger 

E ( sup \\xMHu.)XiiH^)Xy\\l] < Ce-^l^-^l (50) 

\\9\<1 J 

for some rj > and intervals / in the localized regime. 

It would be interesting to find out if the FMM applies to the random displacement model 
under the assumptions of Theorem 9.1. While the FMM uses Lifshitz tail bounds in a similar 
way as MSA, one would need to replace the Wegner estimate by another technical tool, the 
fractional moment bound 

sup E(||x..(if. -E- ter'xyW') < Ce-"!^-^! (51) 

Eel,e>0 

for suitable s G (0, 1). In fact, it would already be a major step to establish finiteness of the 
left hand side of (51) as an a priori bound. We believe that the properties of the RDM which 
went into our proof of a Wegner estimate can also lead to a proof of this a priori bound. 

11.4. Non-generic single site potentials. We have argued above that alternative (i) of 
Theorem 3.1 is the generic case and we have proven our localization results under this as- 
sumption. Nevertheless, we also observed that via formula (9) one gets a rich reservoir of 
examples in which the ground state of the single-site Neumann operator H^ (a) is constant 
outside the support of the potential and thus Eo{a) identically vanishing. When choosing a 
single-site potential of this type in defining the random displacement model, one sees that 
Tsima{H^) = min a{H^i^) = for all choices of u and L. Thus we are not in the fiuctuation 
boundary regime which we exploited above to get the Lifshitz tail bound and thus can not 
prove low energy localization with our methods. 

The determination of the spectral type of Hi^ near zero for this case remains an open 
problem. The fact alone that zero is the deterministic ground state energy for finite volume 
restrictions of H^ does not exclude the possibility of Lifshitz tails, as a single non-degenerate 
eigenvalue does not affect the IDS in the infinite volume limit. On the other hand, the 
ground state eigenfunction in this model is an essentially uniformly spread out extended state 
(differing from a constant only by local random fluctuations). While we don't necessarily 
believe that this could indicate the existence of absolutely continuous spectrum near zero, 
it might lead to non-trivial transport, similar to the existence of critical energies in certain 
one-dimensional random operators such as the dimer model, e.g. [25]. Finding the answer to 
this will require a much better understanding of excited states of the RDM, a task which we 
have managed to systematically avoid in all the results presented above. 



30 klopp, loss, nakamura, and stolz 

References 

[1] M. Aizenman, A. Elgart, S. Naboko, J. Schenker and G. Stolz, Moment analysis for localization in random 

Schrodinger operators, Invent. Math. 163, 343-413 (2006) 
[2] M. Aizenman, F. Germinet, A. Klein and S. Warzel, On Bernoulli decompositions for random variables, 

concentration bounds, and spectral localization, Probab. Theory Related Fields 143, 219-238 (2009) 
[3] M. Aizenman and S. Molchanov, Localization at large disorder and at extreme energies: a elementary 

derivation, Comm. Math. Phys. 157, 245-278 (1993) 
[4] J. Baker, M. Loss, and G. Stolz, Minimizing the ground state energy of an electron in a randomly deformed 

lattice, Comm. Math. Phys. 283, 397-415 (2008) 
[5] J. Baker, M. Loss, and G. Stolz, Low energy properties of the random displacement models, J. Funct. 

Anal. 256, 2725-2740 (2009) 
[6] J. Bourgain and C. Kenig, On localization in the continuous Anderson-Bernoulli model in higher dimen- 
sion. Invent. Math. 161, 389-426 (2005) 
[7] A. Boutet de Monvel, S. Naboko, P. StoUmann and G. Stolz, Localization near fluctuation boundaries via 

fractional moments and applications, J. Anal. Math. 100, 83-116 (2006) 
[8] D. Buschmann and G. Stolz, Two-parameter spectral averaging and localization properties for non- 
monotonic random Schrodinger operators. Trans. Am. Math. Soc. 353, 635-653 (2000) 
[9] R. Carmona, A. Klein and F. Martinelli, Anderson localization for Bernoulli and other singular potentials, 

Comm. Math. Phys. 108, 41-66 (1987) 
[10] R. Carmona and J. Lacroix, Spectral theory of random Schrodinger operators, Birkhauser, Boston, MA, 

1990. 
[11] J. M. Combes, P. D. Hislop and F. Klopp, An optimal Wegner estimate and its application to the global 

continuity of the integrated density of states for random Schrdinger operators, Duke Math. J. 140, 469-498, 

(2007) 
[12] J. M. Combes, P. D. Hislop and S. Nakamura, The i^-theory of the spectral shift function, the Wegner 

estimate, and the integrated density of states for some random operators, Comm. Math. Phys. 218, 

113-130 (2001) 
[13] D. Damanik, R. Sims and G. Stolz, Localization for one-dimensional continuum, Bernoulli-Anderson 

models, Duke Math. J. 114, 59-100 (2002) 
[14] A. Elgart, H. Kriiger, M. Tautenhahn, I. Veselic, Discrete Schrodinger operators with random alloy-type 

potential, http://arxiv.org/abs/1107.2800, 2010 
[15] L. Erdos and D. Easier, Wegner estimate and Anderson localization for random magnetic fields, http: 

//arxiv.org/abs/1012.5185, 2010 
[16] L. C. Evans, Partial Differential Equations, Second Edition, Graduate Studies in Mathematics 19, Amer. 

Math. Soc, Providence, RI, 1998. 
[17] J. Frohlich and T. Spencer, Absence of diffusion in the Anderson tight binding model for large disorder 

or low energy, Comm. Math. Phys. 88, 151-184 (1983) 
[18] F. Germinet, P. D. Hislop, A. Klein, Localization for Schrodinger operators with Poisson random potential, 

J. Eur. Math. Soc. 9, 577-607 (2007) 
[19] F. Germinet, P. D. Hislop, A. Klein, Localization at low energies for attactive Poisson random Schrodinger 

operators, in: Probability and Mathematical Physics, pp. 85-95, CRM Proc. Lecture Notes 42, Amer. 

Math. Soc, Providence, RI, 2007 
[20] F. Germinet and A. Klein, Bootstrap multiscale analysis and localization in random media, Comm. Math. 

Phys. 222, 415-448 (2001) 
[21] F. Germinet and A. Klein, A comprehensive proof of localization for continuous Anderson models with 

singular potentials, http://arxiv.org/abs/1105.0213, 2011 
[22] F. Ghribi and F. Klopp, Localization for the random displacement model at weak disorder, Ann. Henri 

Poincare 11, 127-149 (2010). 
[23] E. M. Harrell, P. Kroger and K. Kurata, On the placement of an obstacle or a well so as to optimize the 

fundamental eigenvalue, SIAM J. Math. Anal. 33, 240-259 (2001) 
[24] P. D. Hislop and F. Klopp, The integrated density of states for some random operators with nonsign 

definite potentials, J. Funct. Anal. 195, 12-47 (2002) 



THE RANDOM DISPLACEMENT MODEL 31 

[25; 

[26; 



[27; 
[28; 
[29; 

[30 

[31 
[32; 
[33; 
[34; 
[35; 

[36 

[37; 

[38; 

[39 

[4o; 

[41 



S. Jitomirskaya, H. Schulz-Baldes and G. Stolz, Delocalization in random polymer models, Comm. Math. 
Phys. 233, 27-48 (2003) 

W. Kirsch, Random Schrodinger Operators - A Course, pp. 264-371 in: Schrodinger Operators, Lecture 
Notes in Physics, Vol. 345, Proceedings Sonderborg 1988, Eds. H. Holden and A. Jensen, Springer, 1989 
W. Kirsch, Wegner estimates and Anderson localization for alloy-type potentials, Math. Z. 221, 507-512 
(1996 ) 

W. Kirsch, An invitation to random Schrodinger operators, pp. 1-119 in: Panor. Syntheses 25, Random 
Schrodinger operators, Soc. Math. France, Paris, 2008 

A. Klein, Multiscale analysis and localization of random operators. Random Schroinger operators, 121- 
159, Panor. Syntheses, 25, Soc. Math. France, Paris, 2008 

F. Klopp, M. Loss, S. Nakamura and G. Stolz, Localization for the random displacement model, http: 
//arxiv.org/abs/1007.2483, 2010. 

F. Klopp, Localization for semiclassical continuous random Schrodinger operators. II. The random dis- 
placement model, Helv. Phys. Acta 66, 810-841 (1993) 

F. Klopp, Localization for some continuous random Schrodinger operators, Comm. Math. Phys. 167, 
553-570 (1995) 

F. Klopp. A short introduction to the theory of random Schrodinger operators. Vietnam J. Math., 
26(3): 189-216, 1998. 

F. Klopp, S. Nakamura, Spectral extrema and Lifshitz tails for non monotonous alloy type models. Comm. 
Math. Phys. 287, 1133-1143 (2009) 

F. Klopp and S. Nakamura, Lifshitz tails for generalized alloy type random Schrodinger operators. Anal. 
PDE 3, 409-426 (2010) 

L. Pastur and A. Figotin, Spectra of Random and Almost-Periodic Operators, Springer, 1992 
M. Reed and B. Simon, Methods of Modern Mathematical Physics, Vol. 4, Analysis of Operators, Aca- 
demic Press, New York, 1978 

P. StoUmann, Caught by Disorder: Bound States in Random Media, Progress in Mathematical Physics 
20, Birkhauscr, Boston, MA, 2001 

G. Stolz, Localization for random Schrodinger operators with Poisson potential, Ann. Inst. H. Poincare 
Phys. Theor. 63, 297-314 (1995) 

I. Veselic, Wegner estimate and the density of states for some indefinite alloy type Schrodinger operators, 
Lett. Math. Phys. 59, 199-214 (2002) 
F. Wegner, Bounds on the density of states in disordered systems, Zeitschrift fiir Physik B 44, 9-15 (1981) 



1 LAGA, U.M.R. 7539 C.N.R.S, Institut Galilee, Universite de Paris-Nord, 99 Avenue J.-B. 
Clement, F-93430 Villetaneuse, France et Institut Universitaire de France 
E-mail address: klopp@math.univ-parisl3.fr 

^ Georgia Institute of Technology, School of Mathematics, Atlanta, Georgia 30332-0160, 
loss@math.gatech.edu 

E-mail address: loss@iiiath.gatech.edu 

^ Graduate School of Mathematical Sciences, University of Tokyo, 3-8-1 Komaba, Meguro- 
Ku, ToKYA, Japan 153-8914, shu@ms.u-tokyo.ac.jp 
E-mail address: shu@ms . u-tokyo .ac.jp 

"^ University of Alabama at Birmingham, Department of Mathematics, Birmingham, Al- 
abama 35294-1170, STOLZ@UAB.EDU 
E-mail address: stolz@uab.edu 



